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Orthonormal polynomials are employed to compute the transfinite diameter of a set 
consisting of two collinear line segments. With polynomials of degree 10, two-place accuracy 
has been achieved. 



The concept of the transfinite diameter t(E) of a 
closed bounded point set E was introduced by M. 
Fekete [5], 1 and the relation between this domain 
functional and conformal mapping and potential 
theory has been stressed by both Fekete and Szogo 
[7]. According to Fekete 's original definition, 



where 



G),- 

r(E)=lim K2 \IV n , 



F»=max|V(2 1 ,2 2 , . . .,g„)|, 

zttE 



(1) 

(2) 



and where V is the Vandermonde determinant. 

The transfinite diameter is known explicitly for a 
number of elementary geometrical configurations, 
but in general, its numerical evaluation is attended 
by considerable difficulty. In a previous paper [4], 
Davis and Rabinowitz showed that complex orthog- 
onal polynomials may be employed profitably for 
its numerical determination in the case of certain 
simply connected domains, a square and nonconvex 
"bean-shaped" domain having been tested. Re- 
cently, Professor Fekete called the attention of the 
author to the possibility of using an analogous 
relationship in the case where E consists of a finite 
number of rectifiable Jordan arcs. Making use of 
very general theorems recently established by 
Fekete and Walsh [6, p. 61], we ma} r assert: 



Let p ri (z)=a n z n -\-. . . , a w >0 (n=0, . . . ) be complex 
polynomials that are orthonormal over E in the sense 
that 



Then we have 



L 



Pm(z)Pn(S)ds = d rr . 



lim (l/a n y= T (E). 



(3) 

(4) 



For two collinear line segments of equal length 
placed symmetrically with respect to 0, say, E: 
-l<x<—a,a<x<l, 0<a<l, the value of r(E) is 
known theoretically and is simply 



T (E)=iVT^. 



(5) 



Figures in brackets indicate the literature references at the end of this paper. 



For two collinear line segments of unequal length, 
its value has been obtained by N. Achiezer (1) and 
can be expressed as the ratio of certain elliptic 
functions. For more than two line segments a 
closed-form value of r{E) is not known to the author. 
The relationship (4) was tested numerically on 
SEAC to see what could be achieved by way of 
accuracy, using single precision codes. In these 
computations, the value a=l/2 was selected, leading 
in (5) to 

t=JV3~0.4330127. (6) 

The computation was carried out by making use of 
a single precision (8 significant decimals) floating 
point orthonormalizing routine developed by J. 
Bram. This multiple-purpose code has already 
been described in [3], The inner products (3) were 
computed by means of a 10-point Gaussian quadra- 
ture rule on each of the two segments (— 1,— |) and 
(|,1). The machine was programmed to print out 
the coefficients of the orthonormal polynomials, as 
well as the values of these polynomials at the Gaus- 
sian abscissas employed. In this way, it was 
possible to monitor the obvious global properties of 
the orthonormal polynomials, as well as to see where 
the accumulated round-off began to vitiate the 
computations. One way in which this was done was 
as follows: The polynomials p n (x), orthonormal over 
the set E\{ — 1,— -J), (1^,1), are alternately even and 
odd. As n increases, the theoretical zero values for 
the alternate coefficients of p n become contaminated 
by round-off, and these "zeros" begin to assume the 
proportions of the nonzero coefficients and of the 
values of the orthonormal polynomials themselves 
(see table 1). It was found that all significance was 
lost when attempt was made to go beyond n=10. 
The last value of (l/a n ) 1/n gave the value of r{E) 
correct to within 0.009. Additional accuracy is 
obtainable from (4) only by employing double- 
precision coding and going beyond w=10. 

Although closed-form expressions for the ortho- 
normal polynomials over E:(— 1, — §),(§, 1) are not 
available, such expressions are available for the 
Tschebyscheff polynomials for E [2, p. 287]. Here 
the even and odd polynomials have a totally different 
structure. Using the latter polynomials as a guide 
(in the theory of domain polynomials these two sets 
frequently behave alike), we can confirm the slightly 
higher values for (l/a n ) 1/n registered in table 2 for 
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Table 2. Computation of transfinite diameter from leading 
coefficients of orthonormal polynomials 



n 


On 


i/y«n 




V0"/«n+2 


1 
2 
3 
4 
5 

6 
7 
8 
9 
10 


1. 3093073 
4.6017900 
6. 1932355 
23.822212 
31. 868611 

125.35317 
167. 38858 
663. 46546 
885. 36970 
3521. 1597 


0. 76376 
. 46616 
. 54454 
. 45264 
. 50041 

. 44700 
. 48120 
. 44389 
. 47048 
. 44191 


0. 45979 
. 43951 
. 44084 
. 43594 
. 43633 

. 43467 
. 43481 
. 43408 






Theoretical value: 
0.4330127 



odd n. We can also conclude from the form of the 
Tschebyscheff polynomials that the ratio (a»/a» +2 ) 1/2 
would be a good estimator for r(E). Table 2 also 
presents these values, and it will be seen that the 
last entry, (a 8 /ai ) U2 , yields r(E) correctly to within 
0.001. 



The author thanks Morris Newman and Patricia 
Ferrant for their assistance with the numerical 
portions of this paper. 
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